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The  High  Altitude  Effects  Simulation  (HAES)  Program  sponsored  by  the 
Defense  Nuclear  Agency  since  the  early  1970  time  period,  comprises  several 
groupings  of  separate,  but  interrelated  technical  activities,  e.g. , ICECAP 
(Infrared  Chemistry  Experiments -Coordinated  Auroral  Program).  Each  of  the 
latter  have  the  common  objective  of  providing  information  ascertained  as  essen- 
tial for  the  development  and  validation  of  predictive  computer  codes  designed 
for  use  with  high  priority  DoD  radar,  communications,  and  optical  defensive 
system. 

Since  the  inception  of  the  HAES  Program,  significant  achievements  and  results 
have  been  described  in  reports  published  by  DNA,  participating  service  labora- 
tories, and  supportive  organizations.  In  order  to  provide  greater  visibility 
for  such  information  and  enhance  its  timely  applications,  significant  reports 
published  since  early  calendar  197^  shall  be  identified  with  an  assigned  MES 
serial  number  and  the  appropriate  activity  acronym  (e.g.,  ICECAP)  as  part  of  the 
report  title.  A complete  and  current  bibliography  of  all  HAES  reports  issued 
prior  to  and  subsequent  to  HAES  Report  No.  1,  dated  5 February  197!|  entitled , 
"Rocket  Launch  of  an  SWIR  Spectrometer  into  an  Aurora  (ICE  CAP  72),"  AFCRL 
Environmental  Research  Paper  No.  ;i 66,  is  maintained  and  available  on  request 
frnm  TTASTAC,  DoD  Nuclear  Information  and  Analysis  Center,  816  State  Street, 

Santa  Barbara,  California  93102,  Telephone  (805)  9&5-0551. 

This  report,  which  is  the  final  report  under  DNA  contract  001-7^1  -C-Olt 6 
is  the  fourteenth  report  in  the  HAES  aeries  and  covers  technical  activities 
perfermed  during  the  period  February  197^  through  February  1975.  Portionc  of 
this  program  were  reported  previously  in  DNA  Report  329 entitled  "ICECAP 
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Analysis:  Mechanisms  for  Energy  Deposit  in  the  Auroral  Ionosphere;"  in  DNA 

i 

Report  3566F  entitled  "ICECAP  Analysis:  Energy  Deposit  and  Transport  In  the 
Auroral  Ionosphere;"  and  in  other  reports  listed  in  the  KAES  bibliography 
and  pre-HAES  bibliography.  The  purpose  of  the  work  reported  herein  was  to 
investigate  the  importance  of  various  high-altitude  energy  input  and  transport 
mechanisms  of  interest  to  ICECAP  experiments;  to  provide  empirical  quantitative 
models  of  these  processes  as  the  data  allow;  to  assist  in  ICECAP  data  inter- 
pretation; and  to  assist  in  planning  and  execution  of  the  ICECAP  program. 

We  extend  our  sincere  thanks  to  Dr.  T.W.  Watt  of  the  Stanford  Research 
Institute  for  supplying  the  Chatanika  radar  data  and  to  Mr.  D.L.  Carr,  Dr. 

R.D.  Sharp  and  Dr.  E.G.  Shelley  of  the  Lockheed  Palo  Alto  Research  Laboratory 
for  providing  the  incident  electron  and  proton  data  obtained  in  the  coordinated 
Auroral-Zone  experiment.  These  data  enabled  us  to  test  the  AURORA  code. 

We  also  express  our  sincere  gratitude  to  Dr.  P.M.  Banks  of  the  University  , 

of  California  at  San  Diego  for  his  assistance  in  the  modification  of  the  AURORA 
code  and  to  Dr.  C.A.  Blank  of  the  Defense  Nuclear  Agency  for  his  technical 
support  and  coordination  of  the  research  effort. 
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A high-altitude  nuclear  explosion  (HANE)  greatly  disturbs  the  ionosphere- 
magnetosphere  system,  enhancing  the  excitation,  ionization  ^nd  temperature  of  the 
upper  atmosphere,  producing  new  chemical  species,  causing  plasma  irregularities, 
and  storing  in  the  magnetosphere  a long-term  source  of  highly- ionizing  radiation. 
The  attendant  emissions  in  the  infrared,  visible,  and  ultraviolet  interfere  with 
mi ssile-t racking  devices.  The  enhanced  ionization,  especially  in  the  irregulari- 
ties, affects  camnuni cat ions,  radar,  and  navigation  systems.  It  is  therefore 
important  to  predict  the  conditions  that  would  ensue  in  the  event  of  a HANE  and 
to  determine  the  effects  on  operational  systems. 

Cue  of  the  most  important  considerations  is  the  distribution  of  the  ionization  in 
the  atmosphere  produced  by  the  relativistic  electrons  released  by  a nuclear  deto- 
nation. The  AURORA  code  has  been  developed  at  IMSC  to  treat  this  problem.  The 
code  includes  the  effects  of  collisions  and  the  converging  magnetic  field.  The 
computations  have  been  verified  by  laboratory  experiments  and  by  geophysical  obser- 
vations at  relatively  low  L values  where  collision  effects  in  the  atmosphere  are 
clearly  dominant  processes.  However,  in  the  auroral  zone,  as  in  the  case  of  a 
HANE,  additional  processes  such  as  plasma  instabilities  ani  electric  fields  may 
affect  the  motions  of  the  precipitating  electrons.  It  is  therefore  important  to 
determine  the  importance  of  these  processes.  Accordingly,  several  tests  of  this 
code  have  been  made  utilizing  the  '•esults  of  coordinated  experiments  wherein  the 
precipitation  electrons  were  observed  with  a low-altitude  satellite  and  the  dis- 
tribution of  the  electron  number  density  in  the  atmosphere,  resulting,  from  the 
incident  electrons,  was  measured  with  the  Chatanika  incoherent  scatter  radrr. 
Results  of  one  test  were  presented  in  a previous  report  (Ref.  1).  New  results  on 
a second  case  are  described  in  Section  2.  in  every  case  the  ionization  profile 
computed  by  the  code  was  in  agreement  with  the  experimental  data,  within  the  accur- 
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acies  of  the  pertinent  parameters.  However,  the  experiments  were  not  sufficiently 
precise  to  provide  a definitive  test  of  the  code. 

Several  independently-developed  codes  are  now  in  use  on  DNA  programs  to  compute  at- 
mospheric effects  of  precipitating  electrons.  A meeting  was  therefore  organized  t 
discuss  the  Clements  of  the  codes  and  to  compare  the  results  computed  by  the  codes 
for  the  same  test  problems.  The  data  presented  at  the  meeting  are  discussed  in 
Section  3. 


Seme  new  applications  cf  the  AURORA  code  involving  the  precipitation  of  fission- 
decay  electrons  in  the  atmosphere  are  also  discussed  in  Section  3. 


; 


Recently,  the  original  version  of  the  VCSD  code  (see  Section  U)  was  incorporated 
in  the  AURORA  code  to  compute  electron  fluxes  in  the  energy  range  1 - 500  eV 
(primaries  and  secondaries)  that  result  from  the  interactions  f precipitating 
electrons  with  tilt  atmospheric  constituents.  This  low-energy  portion  cf  the  code 
has  been  revised,  as  discussed  in  Section  h,  and  can  now  be  reliably  applied  to 
problems  involving  various  chemical-reaction  rates,  atomic-excitation  rates,  and 


optical-emission  rates. 


Finally,  in  Section  5 a wav«- particle  instability  is  discussed  which  may  seriously 
affect  the  distribution  of  the  trapped  electrons  ensuing  from  a nuclear  detonation. 
Among  other  effects,  the  instability  could  -a  ;r  rap::  pitch-angle  diffusion  nt  the 
electrons  at  relatively  low  altitudes,  and  thereby  modify  the  ionization  profiles 
predicted  by  the  existing  computer  codes. 
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2.  ANALYSIS  OF  THE  COORDINATED  EXPERIMENT  IN  AURORAL  ZONE 

in  a previous  study  (Ref.  1),  we  utilized  the  existing  data  from  a coordinated 
radar- satellite  experiment  (January  27,  1972)  to  test  whether  the  AURORA  code 
includes  all  important  theoretical  processes  required  to  calculate  correct  auroral 
electron  energy  deposition  rates  in  the  atmosphere.  The  experiment  involved  the 
DNA  incoherent- scatter  radar  at  Chatanika,  Alaska  (Ref.  2)  and  the  polar-orbiting 
satellite  1971-085A  (Ref.  3). 

The  results  obtained  in  that  investigation  indicated  that  the  theory  incorporated 
in  AURORA  is  adequate  to  account  for  auroral  electron  energy  deposition  rates  in 
the  atmosphere.  However,  due  to  the  large  experimental  uncertainties  and  the  poor 
radar- satellite  geometry  associated  with  that  case,  it  was  recognized  that  addi- 
tional cases  must  be  investigated  in  order  to  obtain  a more  conclusive  result. 

We  have  recently  investigated  another  case  employing  data  from  a radar- satellite 
experiment  conducted  near  1331  UT  on  8 December  1971.  The  same  analytic  procedure 
described  in  Reference  1 was  used  to  analyze  these  data. 

Only  the  high-energy  (>  500  eV)  part  of  the  AURORA  code  was  used  in  this  investi- 
gation. The  effects  due  to  the  modification  of  the  low-energy  part  of  the  code 
are  discussed  in  Section  4.  The  contribution  to  electron  energy  deposition  rates 
below  500  eV  were  estimated  by  extrapolation  of  the  AURORA  results  into  this  region. 
The  correction  amounted  to  only  ~ 5 percent  in  the  region  of  the  peak,  and  ~ 10 
percent  at  higher  altitudes. 

It  will  be  recalled  that  the  high-energy  code  numerically  solves  the  appropriate 
Fokker- Planck,  steady-state  diffusion  equation  to  determine  the  flux  of  electrons 
produced  in  the  atmosphere  by  a given  spectrum  of  incident  auroral  electrons.  The 
solution  takes  into  account  atmospheric  scattering  electron  energy  loss,  and  the 
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mirroring  effect  of  the  geomagnetic  field.  The  rate  at  which  energy  is  deposited 
in  the  atmosphere  by  electrons  is  also  calculated  by  the  program. 

in  the  coordinate!  radar-satellite  experiment,  particle  detectors  on  the  satellite 
measured  the  incident  flux  of  auroral  electrons  at  the  satellite  altitude.  Simul- 
taneously, radar  measurements  were  made  to  determine  the  electron  number  density 
altitude  profile  produced  in  the  E-region  by  the  incident  electrons  interacting 
with  atmospheric  particles.  To  test  the  accuracy  of  the  program,  this  experimental 
E-region  electron  density  profile  is  compared  with  the  theoretically  predicted 
profile  derived  from  electron  energy  deposition  rates  calculated  by  AURORA  for 
the  measured  incident  auroral  electron  flux.  The  theoretical  profile  is  detemined 
from  the  equation  N = /q fa  that  expresses  the  balance  between  electron  production 
and  loss  in  the  E-region.  N is  the  electron  number  density,  a is  the  effective 
recombination  coefficient,  and  q is  the  ion-pair  production  rate  produced  by  the 
auroral  electrons.  All  three  quantities  are  a function  of  altitude,  Available 
experimental  values  of  a pertaining  to  the  auroral  zone  are  used  in  the  equation 
for  N,  and  are  given  in  Ref.  1.  They  are  also  shown  in  Figure  1.  Values  of  q are 

O 

determined  from  energy  deposition  rates  (ev/cm/sec)  obtained  from  AURORA  after 
division  by  35  ev/ion  pair. 

Figure  2 shows  a plan  view  of  the  geometry  between  the  radar  beam  (azimuth  = 279°, 
elevation  = 68  ) and  the  satellite  trajectory  at  795  Km  altitude  mapped  down  geo- 
magnetic field  lines  to  110  Km  altitude,  the  altitude  region  of  major  interest. 

The  radar  beam  and  the  field  line  through  the  satellite  position  at  1331:09  OT 
intersect  the  indicated  altitudes  at  the  locations  shown.  Precipitating  auroral 
electrons  tend  to  J.ie  along  lines  of  constant  geomagnetic  latitude.  The  line  cf 
constant  geomagnetic  latitude  passing  through  the  radar  region  of  measurement  at 
110  km  altitude  would  intersect  the  satellite  trace  approximately  halfway  between 
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Experimental  values  of  the  effective  recombination  coefficient  G.  Omholt’ s values  were  revised 
by  Pcpcff  et  al.  (Ref.  5)  using  a more  recent  value  of  k=  12.5  for  the  ratio  of  total  electron 
production  tc  the  total  iJg  first  negative  band  system.  The  revised  values  are  shown.  Johnson 
and  Bradbury  revised  their  preliminary  va.1  ..es  of  a (Ref.  7)  which  are  published  in  Meyerott  and 
Evans  (Ref. 6).  Their  revised  values,  which  are  reported  by  Popoff  et  al.  (Ref. 5), are  plotted 
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the  8-see  and  9-sec  marks.  Figure  2 indicates  that  at  E- region  heights  the  region 
of  radar  measurements  is  separated  by  distances  of  20-1*0  km  from  the  region  of 
precipitating  electrons  beneath  the  satellite.  Thus,  in  order  to  compare  results 
in  these  two  regions,  as  we  do,  we  must  assume  that  the  characteristics  of  the 
precipitating  electrons  were  approximately  uniform  over  distances  of  that  magnitude. 

Experimental  electron  density  profiles  obtained  with  the  radar  near  the  time  of 
interest  are  shown  in  Figure  3.  Results  are  given  for  two  integration  periods  of 
one  minute  each  and  for  one  integration  period  of  20  sec.  All  three  curves  are 
similar  and  indicate  the  electron  density  profile  in  the  region  of  radar  measurements 
was  fairly  stable  during  the  satellite  pass  near  the  radar  site. 

To  determine  a representative  electron  flux  to  apply  to  the  analysis,  the  particle 
data  were  averaged  over  a four-second  interval  centered  at  1331:09  UT.  Figure  4 
shows  three  incident  electron  flux  energy  spectra  measured  by  differential  energy 
detectors  on  the  satellite  at  1331:09  UT  and  two  seconds  before  and  after  this  time. 
Each  curve  was  derived  from  two  sets  of  detectors.  One  set  was  pointing  toward  the 
zenith  and  the  other  set  was  pointing  95°  away  from  the  zenith.  The  measured  fluxes 
from  the  two  sets  differed  by  factors  of  ^ 2 or  less. 

The  three  curves  in  Figure  4 were  averaged  and  the  resulting  curve  was  used  to  rep- 
resent. the  incident  electron  flux  during  the  experiment.  The  AURORA  code  was  then 
applied  to  calculate  electron  energy  deposition  rates  for  a spectrum  of  this  type, 
assumed  to  be  isotropically  incident. 

Figure  5 shows  the  calculated  electron  density  profile  resulting  from  the  above 
primary  spectrum.  Also  shown  for  comparison  is  the  experimental  profile  for  the 
time  period  1331:00  to  1331:20  OT.  A correction  has  been  applied  to  this  curve 
to  account  for  a small  contribution  due  to  an  auroral  proton  component  which  was 
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ELECTRON  DENSITY  (105  EL/CM3) 

Figure  3-  Experimental  electron  density  profiles  obtained 
during  the  coordinated  radar-stellite  experiment. 
Integration  periods  are  indicated. 
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Fifrure  1,  Incident  electron  energy  spectra  determined  from 
particle  measurements  made  on  board  the  satellite 
durinr  the  coordinated  radar- satellite  experiment. 
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Figure  5-  Comparison  of  experimental  and  calculated  electron 

density  profiles  for  the  coordinated  radar~satellite 
experiment.  The  experimental  profile  has  been  corrected 
for  a small  contribution  from  an  auroral  prot  n component. 
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* present  during  the  experiment  and  detected  by  proton  detectors  on  the  satellite. 

The  correction  amounted  to  less  than  10  percent  at  all  altitudes.  The  procedure 

€ 

for  estimating  the  proton  correction  is  described  in  Reference  1. 

Note  that  the  two  curves  are  in  good  agreement  along  the  low-altitude  edge  of  the 
ionization.  The  peak  values  are  displaced  in  altitude  by  several  kilometers.  At 
the  higher  altitudes  the  calculated  density  values  exceed  the  experimental  values. 
The  discrepancies  noted  between  the  curves  are  large,  yet  are  within  the  limits 
of  uncertainty  introduced  by  the  inaccuracies  in  particle  measurements  and  the 
poor  radar- satellite  geometry. 

Because  of  the  large  displacement  in  altitude  between  the  peaks  of  the  two  curves, 
the  comparison  is  less  satisfactory  for  this  case  than  it  was  for  the  first  case 
where  this  particular  discrepancy  was  small.  A large  fluctuating  latitude  gradient 

S 

i 

was  observed  in  the  present  satellite  data,  but  not  in  the  previous  measurements  of 

* January  27.  Thus  the  averaged  flux  used  in  the  present  analysis  has  a larger 
uncertainty  associated  with  it,  which  may  account  for  the  less  favorable  comparison 
between  the  experimental  and  theoretical  peak  positions.  The  discrepancy  may  also 
be  due  to  an  accelerating  electric  field  acting  on  the  incident  electrons  at 

; 

higher  altitudes.  Such  a field  would  cause  a downward  displacement  of  the  observed 

peak  position  as  discussed  in  Ref.  1. 

i 

1 

t 

We  have  now  analyzed  a total  of  two  cases,  in  both  instances,  the  theoretical  and 

J 

experimental  results  are  consistent  with  each  other,  within  the  uncertainties 
involved  in  the  analysis.  Because  these  uncertainties  are  large,  additional  cases 

* should  be  investigated  to  provide  a set  of  results  which  will  lead  to  a more 
definite  conclusion  about  the  adequacy  of  the  theory  incorporated  in  the  AURORA 

* 

program. 


17 


Section  3 


APPLICATIONS  OF  AURORA  CODE 
3.1  COMPARISON  OF  ENERGY  DEPOSITION  CODES 

On  June  1*+,  197*+,  a meeting  was  held  at  Lockheed  Missiles  and  Space  Co.,  Inc. 

(IKSC)  to  compare  the  results  of  various  computer  programs  which  calculate  electron 
transport  in  the  atmor^.-iere.  In  attendance  were  Major  Larry  Doan  and  Dr.  Charles 
Blank,  DNA;  Dr.  Peter  Banks,  University  of  California,  San  Diego;  Dr.  Andrew  Nagy, 
University  of  Michigan;  Dr.  Timothy  Coffey  and  Dr.  D.  J.  Strickland,  Naval  Research 
Laboratory;  Ellis  Hyman,  Science  Applications,  Inc.,  McLean,  Virginia;  Nicholas 
Winter,  Aerospace  Corporation;  Dr.  0.  Manley  and  Dr.  Henry  Smith,  Visidyne,  Inc.; 

Dr.  T.  L.  Stephens,  General  Electrlc-TEMPO ; Dr.  Douglas  Archer  and  Dr.  Paul  Tarr, 
Mission  Research  Corporation;  and  Dr.  Martin  Walt,  Lockheed  Palo  Alto  Research 
Laboratory. 

prior  to  the  meeting  the  input  electron  fluxes  for  two  sample  problems  were  dis- 
tributed to  various  groups  and  the  individuals  were  instructed  to  calculate  the 
total  energy  deposit  as  a function  of  altitude,  the  spectrum  of  the  upward  and 
downward  moving  electrons  at  several  altitudes,  the  energy  production  of  various 
excited  states  in  atmospheric  constituents,  and  concentrations  of  selected  species 
at  various  times.  Calculations  were  done  by  IMSC  using  a combined  IMSC-UCSD  code, 
Mission  Research  Corporation  (MRC)  using  the  AHTIC  program,  Naval  Research  Labor- 
atory (NRL),  and  Visidyne  (VIS).  Each  of  these  calculations  involved  quite  different 
mathematical  methods  and  approximations.  The  NRL  code  has  been  revised  since  the 
meeting  was  held.  The  NRL  results  presented  at  the  meeting  are  therefore  designated 
* NRL(OLD).  Later  results,  which  are  discussed  subsequently,  are  designated  NRL(NEW), 
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Comparisons  of  the  results  for  the  total  energy  deposit  as  a function  of  altitude 
are  shown  in  Figures  6 and  7-  Problem  1 consisted  of  an  isotropic  input  flux  at 
250  Ion  altitude  of  the  form  f(E)  = exp  (-E/lO  kev)  with  a high  energy  cutoff  at 
50  kev.  Problem  2 involved  the  experimental  spectrum  shown  in  Figure  8.  This 
spectrum  of  electrons  was  assumed  to  be  isotropically  incident  at  200  km  altitude 
along  magnetic  field  lines  having  a dip  angle  of  77°.  The  CIRA  1965  atmosphere 
was  used  in  both  problems. 

Although  the  MRC  and  IMSC  results  are  in  fair  agreement  with  each  other,  the  NRL 
and  Visidyne  results  are  substantially  different,  particularly  for  Problem  1.  In 
general  the  Visidyne  curve  shows  less  energy  deposited  over  most  altitudes  while 
the  NRL  results  exhibit  larger  energy  deposition  at  high  altitudes  and  less  pene- 
tration to  extreme  ranges. 

Drs.  Manley  and  Smith  of  Visidyne  stated  that  they  had  inadvertently  used  a low 
energy  cutoff  of  6 kev  in  their  calculations.  This  cutoff  would  substantially 
reduce  the  input  energy  and  could  explain  why  their  energy  deposit  was  less  than 
the  other  calculations. 

A detailed  comparison  between  the  NRL  and  LMSC  flux  values  indicated  that  the 
electrons  scattered  in  pitch  angle  more  rapidly  with  the  NRL  code.  This  con- 
clusion was  most  apparent  from  the  net  upward  and  downward  fluxes  as  calculated 
by  the  two  methods.  Using  the  NRL-OLD  code  in  Problem  1 the  flux  at  10  kev  was 
almost  isotropic  after  penetrating  to  110  km  in  contrast  to  the  IMSC  code  in 
which  the  downward  flux  was  about  ten  times  the  upward  flux.  Drs.  Strickland  and 
Coffey  pointed  out  that  their  method  of  solving  the  transport  equation  becomes 
difficult  at  high  energies  and  that  the  scattering  treatment  might  need  some  modi- 


fication there. 
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TOTAL  ENERGY  DEPOSITION  (eV/CM  -SEC) 


Calculated  total  energy  deposition  as  a function  of  altitude 
for  an  electron  spectrum  of  the  type  F(E)  = exp  (-E/10  key) 
isotropically  Incident  at  250  km  altitude.  Shown  are  curvea 
obtained  by  Lockheed  Miasiles  and  Space  Company,  Inc.  (LMSC) 
Misalon  Research  Corporation  (MRC),  Naval  Reaearch  Laboratory 
(NRL),  Visidine , Inc.  (VIS),  and  Berger,  Seltser  and  Maeda  (BSM) 
(Ref.  n). 


TOTAL  ENERGY  DEPOSITION  (eV/CM^-SEC) 


Figure  7 . Calculated  total  energy  deposition  as  a ♦’unction  of  altitude  for  th" 
experimental  electron  spectrum  shown  in  Figure  8,  assumed  tc  be  isotropically 
incident  at  200  km  altitude  along  magnetic  field  lines  having  a dip  angle  of 
77°.  Shown  are  curves  obtained  by  Lockheed  Missiles  and  Space  Company,  Inc. 
(IMSC),  Mission  Research  Corporation  (MRC),  Naval  Research  Laboratory  (NRL), 
and  Visidyne,  Inc.  (VIS). 
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Since  the  NRL  energy  deposition  and  electron  fluxes  did  not  agree  with  those  of 

MRC  and  IMSC-UCSD,  there  was  no  point  in  comparing  the  NRL  excitation  rates  or 

chemistry  outputs  with  the  other  calculations.  It  was  not  possible  to  compare 

the  upward  and  downward  fluxes  of  the  IMSC-UCSD  calculations  with  MRC  since  the 

MRC  ARTIC  code  did  not  print  out  these  fluxes.  The  IMSC-UCSD  program  normally 

calculates  excitation  rates  of  major  auroral  emissions  (i.e.  6300  X,  5577  X,  Kg+  2N) 

rather  than  the  levels  requested  by  DNA.  However,  for  Problem  2 the  code  was  modi- 

3 + 

fied  slightly  to  give  the  production  rate  of  the  A I u state  of  Ng.  This  result 
was  in  poor  agreement  with  the  MRC  result,  but  the  discrepancy  was  quickly  identified 
as  coming  from  the  use  of  different  sets  of  Green  and  Stolarski  cross-section 
parameters  by  the  MRC  and  IMSC-UCSD  codes. 

As  mentioned  above,  subsequent  to  the  meeting  the  NRL  code  was  modified  and  rerun 
for  Problem  1.  The  new  results  (unpubi;shed)  are  designated  NRL(NEW)  in  Figure  6. 
Also  some  results  by  Berger,  Seltzer  and  Maeda  (Reference  11)  have  appeared  that 
can  be  compared  with  the  problem  1 results.  Though  their  work  was  not  connected 
with  the  DNA- electron  deposition  codes  exercise,  one  of  their  cases  was  identical 
with  the  specifications  for  problem  1.  The  results  are  labelled  BSM  in  Figure  6. 
Berger  et  al.  used  an  isotropic  input  spectrum,  which  was,  however,  normalized  to 
a horizontal  plane.  Their  atmospheric  model  was  the  same  as  specified  in  problem 
1.  Berger  et  al.  used  a Monte  Carlo  method  which  yields  relatively  precise  results 
and  must  take  account  of  scattering  at  all  angles.  (The  Monte  Carlo  method  leads 
to  long,  tedious  computations,  and  is  probably  not  as  suitable  as  the  Fokker- Planck 
method  (as  in  AURORA)  for  DNA's  operational  requirements.)  Comparison  of  the 
results  shows  that  the  new  NRL  calculations  still  predict  slightly  higher  energy- 
deposition  values  at  the  higher  altitudes  although  the  disagreement  is  not  large. 

It  is  gratifying  that  three  independent  calculations,  IMSC,  MRC,  and  BSM  yield 
nearly  identical  results. 
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3.2  ENERGY  DEPOSITION  RESULTING  FROM  HANE 


Two  computer  runs  were  made  with  the  AURORA  computer  code  to  obtain  energy  deposi- 
tion results  for  two  HANE  problems  supplied  by  Dr.  John  Ise  of  General  Research 
Corporation  at  Sants.  Barbara,  California.  The  problems  involved  beta  particles 
produced  in  a high-altitude  nuclear  detonation  and  injected  along  geomagnetic 
field  lines  into  the  atmosphere.  Both  problems  were  identical  except  for  the  dip 
angle  made  by  the  field  liras,  70.1°  (L  = 2.9)  in  one  case  and  32.3°  (L  = 1.1)  in 
the  other  case. 


The  CIRA  atmosphere  was  used  in  the  calculation.  The  minimum  and  maximum  energy 
limits  imposed  on  the  calculation  were  1 keV  and  10  Mev.  The  incident  beta  particle 
spectrum  was  specified  by  the  formula,  I.36  x 10  exp(-E)  (cm  -sec-ster-MeV)  , 
where  E is  in  Mev  units.  This  spectrum  was  assimed  to  be  isotropically  incident 
along  magnetic  field  lines  at  150  km  altitude.  The  incident  flux  crossing  a surface 
perpendicular  to  the  field  lines  may  be  easily  calculated  and  is  found  to  be  4.35  x 
10'  cm  sec  (E  > 1 keV).  The  magnitude  of  the  incident  flux  corresponds  to  a 
realistic  situation,  and  is  an  estimate  of  the  beta  particle  flux  impacting  on  the 
atmosphere  at  150  km  when  a 5-MT  nuclear  device  is  detonated  at  high  altitude. 

The  results  of  the  computer  calculations  are  given  in  Figure  9,  which  shows  the 
energy  deposition  rates  as  a function  of  altitude  for  both  problems. 

Subsequent  to  the  calculations  made  by  IMSC,  results  for  the  same  problem  in  which 
the  magnetic  dip  angle  = 70.1°  were  also  calculated  by  Dr.  Ise,  using  a General 
Elect ric-TIMPO  computational  model  (Ref.  12),  and,  presumably,  the  same  specified 
conditions.  Electron  scattering  is  ignored  in  the  model,  and  the  reflection  of 
electrons  by  the  converging  geomagnetic  field  is  not  taken  into  account.  The 
results  are  presented  as  curve  B in  Figure  10.  The  corresponding  results  obtained 
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Fit^ure  <p.  Fjaerf^-  deposited  as  o function  of  altitude  in  the  atmosphere 
by  beta  particles  from  two  hypothetical  hlrh-altitude  nuclear 
detonations,  one  aasumed  to  occur  where  the  magnetic  dip 
anr,le  = 32.3°  sad  the  other  assumed  to  occur  where  the  dip 
anp;le  - T’.l  . See  text  for  details  of  the  assumed  nuclear 
source. 


ALTITUDE  (KM) 


! 

! 


Figure  10.  Calculated  electron  energy  deposition  as  a function  of  altitude. 

Curves  A and  C were  obtained  by  LhSC  and  curve  B was  obtained 
with  a General  Electric-TJMPO  computational  model  (Ref.  12). 
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previously  for  this  problem  (see  Figure  4)  by  IMSC  are  given  by  Curve  A.  Note 

2 

the  energy  deposition  is  expressed  in  MeV  units  in  Figure  9 and  in  Mo  units  in 
Figure  10. 

A comparison  of  the  two  solutions,  supposedly  obtained  for  the  same  problem, 
indicates  discrepancies  that  are  significant  enough  to  require  an  explanation. 

In  the  previous  IMSC  calculation,  the  given  spectrum  of  electrons  was  assumed  to 
be  isotropic  in  pitch  angle  at  the  incident  altitude  of  150  km.  However,  in  the 
GE-TEMFO  formulation  of  the  problem,  the  pitch  angle  (9)  distribution  of  the  incident 
spectrum  turned  out  to  be  proportional  to  l/cos  9 rather  than  isotropic.  (The  GE- 
TEMFO  angular  distribution  is  appropriate  for  a layer  of  fission  fragments  located 
near  the  top  of  the  atmosphere,  while  the  IMSC  angular  distribution  is  better  for 
fission  fragments  emitting  beta  particles  well  above  the  atmosphere.)  Because  the 
GE-TEMFO  distribution  contains  more  electrons  at  large  pitch  angles,  electrons  do  not 
penetrate  as  far,  on  the  average.  As  a consequence,  electron  energy  is  deposited 
at  higher  altitudes  (curve  B)  as  compared  to  the  results  obtained  by  IMSC  (curve  A). 

To  make  the  comparison  more  direct,  the  same  problem  was  run  again  with  the  IMSC 
AURORA  code,  except  that  the  input  angular  distribution  was  made  proporatlonal  to 
l/cos  9 rather  than  isotropic.  Curve  C in  Figure  10  shows  the  results  of  this 
calculation. 

The  discrepancies  between  C and  the  other  curves  are  attributed  to  the  following. 

For  particles  injected  at  150  km  with  a l/cos  9 distribution,  a significant  fraction 
( ~20  percent)  are  reflected  by  the  magnetic  field  before  they  penetrate  very  far 
into  the  atmosphere;  another  fraction  ( -10  percent)  are  backscattered.  Hence,  the 
total  energy  loss  (given  by  the  integral  under  the  curves)  is  appreciably  less  for 
C than  B.  The  altitude  dependence  of  c is  approximately  the  same  as  for  A because, 
once  the  particles  penetrate  far  enough  into  the  atmosphere,  their  pitch  angles  are 
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isotropized  by  scattering.  Other  differences  may  result  from  the  approximation 
dE/dx  = -2  MeV/gram  cm"2  used  in  the  GE-TEMPO  code,  but  are  believed  to  be  minor 


4.  MODIFICATION  OF  AURORA  CODE 
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The  original  version  of  the  low-energy  (1  - 500  eV)  code  that  has  been  incorporated 
in  the  AURORA  code  is  described  in  detail  by  Banks,  et  al.  (T  f.  13).  This  code 
is  based  upon  solutions  to  continuity  equations  which  take  into  account  transport 
of  electrons  along  magnetic  field  lines  as  well  as  the  effects  of  elastic  and 
inelastic  collisions  with  ambient  neutral  and  Ionized  gases  of  the  upper  atmosphere. 
The  Low  Energy  Auroral  Electron  Code  (LEAURC-l)  provides  values  of  upgoing  and 
downgoing  hemispherical  electron  fluxes  for  energies  1 - 500  eV  over  a preselected 
range  of  altitudes.  Starting  from  a prescribed  incident  electron  flux,  the  high- 
energy  portion  of  AURORA  computes  electron  fluxes  as  a function  of  pitch  angle  and 
altitude  for  electron  energies  greater  than  500  ev.  These  results  are  then  trans- 
mitted to  LEAURC-I  where  the  low  energy  flux  computations  are  made,  in  addition, 
LEAURC-I  carries  out  the  final  calculation  of  atmospheric  auroral  emissions  and 
individual  species  ionization  rates  using  internally  generated  cross-sections  and 
a standard  model  atmosphere. 

Although  LEAURC-I  has  been  applied  successfully  to  practical  problems,  the  code  has 
exhibited  several  undesirable  features.  The  elimination  of  these  features  has 
led  to  the  development  of  an  advanced  version  of  the  code  (LEAURC-II)  which  should 
prove  to  be  substantially  less  costly  to  run  while  giving  more  accurate  results. 

One  of  the  most  important  changes  incorporated  into  LEAURC-II  hss  been  the  elimi- 
nation of  all  nonessential  cross-section  calculations.  As  before,  the  low  energy 
code  obtains  the  data  base  (atmospheric  densities,  electron  fluxes,  altitude 
increments,  etc.)  from  the  high  energy  code.  Following  a short  series  of  internal 
computations  (to  include  secondary  electrons  for  energies  above  500  eV) , the  low 
energy  code  is  prepared  to  begin  its  computations  of  the  upgoing  and  downgoing 
fluxes  in  the  highest  energy  bin  (~  500  eV).  To  make  these  computations,  cross 
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sections  for  elastic  scattering,  inelastic  scattering,  and  scattering  probabilities 
are  read  from  a separate  cross-section  disc-file,  in  addition,  since  discrete  de- 
excitation and  continuum  ionization  processes  are  present,  two  electron  energy 
transfer  matrices  are  provided  to  guide  the  placement  of  cascading  electrons  into 
the  proper  energy  bins.  Finally,  the  various  inelastic  cross  sections  for  0,  0 2 
and  Ng  are  used  to  compute  38  separate  optical  excitation  rates. 

in  the  earlier  version  of  the  low  energy  code  (LEAURC-l),  most  of  the  cross-sections 
were  obtained  either  from  analytical  approximations  to  experimental  data  or  tabu- 
lated data.  Thus,  in  successive  applications  of  the  code,  redundant  computations 
were  invariably  made.  The  newer  veraion  of  the  code  (LEAURC-Il)  avoids  this  redun- 
dancy through  the  use  of  a cross  section  data  file  constructed  to  match  the  demands 
of  the  computational  code.  Although  the  use  of  this  file  introduces  a substantial 

l 

time-lag  between  computations  of  the  electron  flux  in  successive  energy  bins,  such 
a delay  does  not  appear  to  be  a handicap  if  a multi -pro cess or  computer  is  used. 

in  addition  to  its  cross-section  data  file  the  LEAURC-II  code  also  incorporates  an 
improved  energy  mesh,  previously,  an  arbitrary  mesh  was  specified  which  included 
a msnber  of  discontinuities  in  the  mesh  width,  A E,  at  selected  points.  The  size  of 
AE  was  chosen  to  accommodate  a fine  mesh  (~  1 eV)  at  low  energies  (E  < i,/  eV), 
expanding  to  a 20  eV  width  at  500  eV.  The  fluxes  computed  with  this  mesh  are 
regarded  as  being  reasonably  accurate  but  subject  to  systematic  variation  near  the 
transition  energies  where  AE  is  abruptly  changed.  In  LEAURC-II  a new  method  of 
choosing  the  bin  energies  has  been  adopted  so  that  AE  changes  continuously  as  a 
function  of  energy.  Although  more  energy  bins  are  required  (and,  hence,  more 
computations)  a substantial  Improvement  in  the  computed  fluxes  should  result. 


As  a final  note,  during  the  past  year  a number  of  new  data  have  become  available  for 
certain  excitation  and  scattering  cross-sections.  These  newer  data  have  been  incor- 
porated into  the  basic  data  base  and  should  enhance  the  accuracy  of  the  calculations. 
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The  revised  code  was  used  to  compute  the  energy-deposition  profile  in  the  atmos- 
phere resulting  from  the  spectrum  of  precipitating  electrons  shown  in  Figure  8 
(Test  Problem  2) . The  profile  was  found  to  be  indistinguishable  from  the  IMSC 
results  shown  in  Figure  7 at  altitudes  below  about  170  km.  Above  170  km  the 
results  departed  slightly,  with  the  new  profile  indicating  somewhat  higher  energy 
deposition. 
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Section  5 


INTERACTION  OF  FISSION-DECAY  ELECTRONS  WITH  ION-CYCLOTRON  WAVES 

5.1  INTRODUCTION 

The  interaction  of  relativistic  electrons  with  electromagnetic  ion-cyclotron  (ion- 
EMC)  waves  has  been  suggested  as  an  important  trapped- electron  loss  mechanism 
(Ref.  14,  Ref.  15),  and  an  important  trapping  mechanism  for  artificial  rsdistion 
belts  (Ref,  l6) . The  distinction  between  the  two  cases  is  that  the  trapped  elec- 
trons can  be  precipitated  while  absorbing  energy  from  the  ambient  waves,  but  a 
narrow  beam  of  electrons  concentrated  in  the  loss  cone  is  unstable  to  generation  of 
ion-EMC  wavea. 

The  interaction  of  electrons  with  ion-IMC  waves  may  be  important  whenever  an  ab- 
normally high  number  of  relativistic  electrons  is  present  in  the  magnetosphere. 

It  could  be  relevant  to  the  electron  precipitation  problem  and  must  be  examined 
because  it  is  one  mechanism  that  haa  not  been  included  in  the  AURORA  code,  in 
the  following  section  the  general  theory  of  wave-particle  intersetions  in  a finite- 
temperature  plasma  is  discussed.  The  succeeding  section  describes  a particular 
Instability  and  its  relevance  to  the  artificial  electron  belt  problem  and  the 
electron  precipitation  problem. 

5.2  THE  FINITE-TEMPERATURE  RESONANCE  CONDITION 

The  well-known  electromagnetic-mode  resonance  condition  is 

n 

p„  /me  = VY  ' “S*  (1) 

where  p is  the  component  of  momentum  slong  the  magnetic  field  direction,  m is 

the  electrons'  mass,  y is  the  relativistic  "Lorentz"  factor,  and  n ia  the  electron 
' e 
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gyrofrequency.  The  wave  is  characterized  by  the  frequency  uj,  which  is  positive 
for  right-handed  (whistler-mode)  waves  and  negative  for  left-handed  (ion- EMC) 
waves,  and  the  phase  velocity,  V . The  resonance  condition  with  ion-5MC  waves 
is  most  favorable  near  the  ion  gyrofrequency,  where  simple  cold  plasma  theory 
gives  the  dispersion  equation 
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where  is  the  ion  gyrofrequency  and  0^  is  the  plasma  frequency.  For  u>  near  -0^ 
the  phase  velocity  becomes  very  small,  and  the  resonance  condition  (Eq.  1 ) can 
be  satisfied  for  moderate  momenta. 

The  equation  ( 2)  above  is,  however,  not  complete.  Near  the  gyrofrequencies  the 
dispersion  equation  is  substantially  modified  by  the  effects  of  plasma  temperature. 
The  correct  dispersion  equation  is  (Ref.  17,  18,  19  ) 

•V  - v/>  + Mp  [mfr  * Vtr  vJir)  + VI  - 0 <3> 

all  species  11 

where  n and  u)  refer  to  the  gyro-  and  plasma- frequency  of  a particular  constituent 

C P rp 

of  the  plasma.  The  so-called  temperature  anisotropy  is  A_=  (7=^ 1)  , and  F is 

1 il 

defined  as 


F(x) 


(4) 


The  plasma  dispersion  function,  Z(x)  has  been  tabulated  by  Fried  and  Conti  (Ref.  20). 
Near  the  ions*  gyrofrequency  the  "thermal"  part  of  the  distribution  predominates, 
and  the  anisotropy  can  be  assumed  zero. 
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A slightly  simplified  dispersion  equation  can  be  written  in  terms  of  a single 
Maxwellian  tempterature  T: 


<u(c  - V ) + V -~sr-  F ( V 

9 ' 9 L p u)+0  ' u>  9 


) = o 


The  temperature  must  be  related  to  the  "beta"  of  the  plasma  - the  ratio  of  particle 
energy  to  magnetic  field  energy.  We  can  define  three  parameters; 


f)  2 N m c2 


2 

3 B2/8n 


E = ~6<  .000543 
“i 


(The  definition  of  & here  differs  by  a factor  of  ^ from  the  conventional  definition.) 
Substituting  the  resonance  condition  (Eq.  1 ) in  Equation  (5), and  omitting  terms  of 


onier  <u/^e  gives 


Note  that  in  the  limit  T 0,  F approaches  unity  and  Equation  ( 7)  reduces  to  the 
correct  approximation 


P r ne(uj+ni)  |'2  Cne3(a>+ci)3‘ 


But  for  finite  temperatures,  as  u)  approaches  F(x)  is  nearly  The  resonance/ 


dispersion  equation  (7)  then  becomes 

(!if  . (M2 , i 

\mc}  4 n2&  \mc/  y n. 


+ 1 ~ 0 


The  quadratic  solution  is 


i U t/z  „ ri  yi1/2  (10. 

mc  LQ  a.^)  J LQWi  <*>/J 

Equations  (8)  snd  (10)  tell  us  that  p^  becomes  very  large  as  <u  approaches  either  0 

or  ; we  can  deduce  that  there  must  be  a minimum  value  of  p^  for  which  resonance 

is  possible.  That  minimum  should  be  expected  to  occur  near  the  region  where  F(>:) 

2 

switches  from  2x  to  1,  or  near 


^P 


For  small  values  of  0,  Equations  (8)  and  (11)  give 


(1)  ill  M) 

\Q/  min  </Q 

Figure  11  shows  the  minimum  resonant  electron  momenta  (computed  from  Equation  7), 
together  with  the  corresponding  minimum  icn-EMC  wave  phase  velocities  (the  dimension- 
less quantity 


u . 

min  \rac  ymi 


has  been  plotted  there), 


Values  of  Q are  generally  of  the  order  100  inside  the  plssmasphere  and  1 or  less 
outside.  Near  the  plssmapause  0 may  be  about  0.1  during  moderately  disturbed  condi- 
tions; in  that  esse  the  limiting  momentum  is  sbout  1.5  MeV.  A reasonable  value 
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for  6 during  quiet  times  right  be  .001,  which  gives  a limiting  resonant  momentum 
of  about  ,7  MeV.  The  minimum  resonant  momentum  is  not  very  sensitive  to  g.  That 
the  electron  momentum  has  to  be  at  least  .7  MeV  means  that  a large  percentage  of 
fission-beta  electrons  might  participate  in  interactions  with  ion-EMC  waves. 

5.3  THE  AMPLIFICATION  OF  ION-EMC  WAVES  BY  RELATIVISTIC  ELECTRONS 

An  amplification  rate  for  ion-EMC  waves  excited  by  a narrow  stream  of  electrons  is 
readily  derived,  it  is  (Ref.  16) 


where  is  that  part  of  Q contributed  by  electrons  in  the  stream.  For  a stream 

of  electrons  concentrated  within  the  loss  cone  a convenient  distribution  function  is 


(u,,-  M-  " 2UC  + 2u  ) m 3/2  2 

f ( cos  a,  p)  = f(u,j  p)  « 3 (-2Ft-^  exP(-P  /^V1^ 

(1  - ) e 


(1*0 


The  amplification  rate  for  that  distribution  is 


“ ~ '•/"TT  I J-?“  ? 21 
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3/2 


exp  (-p  /2m  T ) 

|1  c c 


(15) 


Equation  (15)  represents  a continuum  spectrum  of  growing  waves,  with  a well  defined 
fastest  growing  mode.  Numerical  solutions  for  the  largest  possible  value  of  Z have 
been  constructed,  using  Equations  ( 7)  and  (15);  the  results  are  plotted  in  Figure  12. 
It  was  nccessai*y  to  distinguish  between  the  temperatures  of  the  ions  and  the  fast 

P 

electron  stream;  in  the  figure  B merely  denotes  the  dimensionless  product  2QT  /m  c 

c c 0 

n 

- to  he  compared  with  B = POT^/m^ . The  effective  temperature  of  fission-beta 
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Figure  12.  Maximum  amplification  rates  for  ion-EMC  waves  excited  by 

relativistic  electrons  for  several  values  of  the  electrons' 
effective  temperature,  denoted  by  2QTrymf.c  . 


electrons  injected  by  a nuclear  explosion  is  sbout  2 MeV,  which  leads  to  a value 
of  about  10^  for  6g  (Q  » 100).  For  most  reasonable  values  of  $ the  maximum  ampli- 
fication rate  is  nearly  Zmax  ^ Slot'  Evidently  the  conditions  are  generally  very 
favorable  to  this  type  of  instability. 

Of  course  Q and  Q„  , vary  along  a field  line,  with  a maximum  expected  near  the 
equator.  But  the  variation  in  is  probably  not  strong,  because  the  density 
of  the  electron  stream  is  approximately  inversely  proportional  to  B,  and  therefore 
compensates  for  the  variation  of  B in  the  denominator  of  Equation  (6a).  We  may 
suppose  that  the  instability  can  take  place  even  at  rather  low  altitudes. 


The  effects  of  the  electron  stream-EMC  wave  instability  are  easy  to  imagine.  We 
have  the  well  loiown  conservation  of  momentum  in  a reference  frame  moving  with  the 
phase  velocity  of  the  wave;  expressed  as  a mathematical  formula: 


p ftp  + (p  - myV  )6p.  * 0 

F±v*i.  ||  9 II 


Together  with  the  resonance  condition  (Eq.  1)  this  gives 


5pj.  ne  Qe 

— i-  ~ ~ ^ -2000 

6^  % 


Apparently  a large  change  in  the  transverse  momentum,  or  pitch-angle,  is  accom- 
plished with  a slight  expenditure  of  energy.  That  the  right  hand  side  is  negative 
confirms  that  a stream  of  electrons  concentrated  in  the  loss  cone  can  give  up 
energy  to  the  waves  while  being  deflected  to  larger  pitch  angles.  Equation  (17) 
alsc.  means  that  a trapped  distribution  of  relativistic  electrons  can  absorb  energy 


from  ion-EMC  waves  and  suffer  decreases  in  their  pitch  angles. 
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The  mechanism  described  here  can  lead  to  enhancement  of  artificial  trapped  electron 
belts,  it  could  also  be  an  important  cause  of  rapid  pitch-angle  diffusion  in  an 
intense  stream  of  precipitating  electrons.  The  diffusion  could  take  place  at  a 
relatively  low  altitude  (above  the  atmosphere),  and  must  be  considered  a potential 
source  of  error  whenever  AURORA  or  other  electron  energy  deposition  codes  are 
applied  to  relativistic  electrons. 
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